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Abstract —This study presents novel methods for computing 
fixed points of positive concave mappings and for characterizing 
the existence of fixed points. These methods are particularly im¬ 
portant in planning and optimization tasks in vrireless networks. 
For example, previous studies have shown that the feasihillty of 
a network design can he quickly evaluated hy computing the 
fixed point of a concave mapping that is constructed based on 
many environmental and network control parameters such as 
the position of base stations, channel conditions, and antenna 
tilts. To address this and more general problems, given a positive 
concave mapping, we show two alternative but equivalent ways 
to construct a matrix that is guaranteed to have spectral radius 
strictly smaller than one if the mapping has a fixed point. This 
matrix is then used to build a new mapping that preserves 
the fixed point of the original positive concave mapping. We 
show that the standard fixed point iterations using the new 
mapping converges faster than the standard iterations applied 
to the original concave mapping. As exemplary applications of 
the proposed methods, we consider the problems of power and 
load estimation in networks based on the orthogonal frequency 
division multiple access (OFDMA) technology. 

I. Introduction 

Problems that can be posed as that of finding fixed points 
of standard interference mappings are ubiquitous in commu¬ 
nication systems m-Ei, and, in particular, in planning and 
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optimization of networks based on the orthogonal frequency 
division multiple access (OFDMA) technology ll-Ha. In 
many of these applications, the mappings are positive concave 
mappings, which are a strict subclass of standard interference 
mappings Q. 

For example, by using standard interference coupling mod¬ 
els that are widely used in the literature ifTTIl . the studies in ||4| 
and El consider a very particular case of a positive concave 
mapping for the problem of load estimation in long-term 
evolution (LTE) networks. The fixed point of that mapping, 
if it exists, indicates the bandwidth required by each base 
station to satisfy the data rate requirements of users. With 
this knowledge, we can evaluate the feasibility of a network 
design by verifying whether the required bandwidth does not 
exceed the available bandwidth. However, especially in large- 
scale planning, computation of the fixed point may require 
time-consuming iterative methods. Therefore, the development 
of fast tools to ensure the existence of a fixed point before 
starting a time-consuming iterative process is of high practical 
relevance to network designers and to algorithms for self¬ 
organizing networks. 

In the above-mentioned load estimation problem, existence 
of a fixed point is fully characterized by the spectral radius 
of a matrix that is easily constructed from the associated 
concave mapping We can also use this matrix to build 
an affine mapping having as its fixed point a vector that gives 
a lower bound of the network load. The main advantage of 
working with affine mappings in finite dimensional spaces is 
that computation of their fixed points reduces to solving simple 
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systems of linear equations, so we may easily obtain in this 
way a certificate that the current network configuration is not 
able to serve the demanded traffic. 

The first objective of this study is to show that, by using the 
concept of recession or asymptotic functions in convex analysis 
HU, m, the technique used in 0 for the construction of 
the above-mentioned matrices (hereafter called lower bounding 
matrices) admits a simple extension to general positive concave 
mappings. This extension has been motivated by recent results 
in power estimation in LTE networks HI, El, which deal with 
mappings different from that considered in a, 0 . Concave 
mappings are also common in many applications in different 
fields 03, so the results of this study are relevant for appli¬ 
cations outside of the wireless domain. We show alternative 
construction methods for lower bounding matrices that are 
very simple in many applications, including those originally 
considered in 0. We also prove that the spectral radius of 
lower bounding matrices of general concave mappings gives 
a necessary condition for the existence of fixed points. For 
some particular concave mappings, this condition is shown to 
be sufficient. 

The second objective of this study is to develop an accelera¬ 
tion method for the standard fixed-point iteration described in 
ifn when applied to concave mappings. More specifically, we 
combine the lower bounding matrix and the original positive 
concave mapping to generate a new mapping that has the 
same fixed point of the original concave mapping. By applying 
the standard fixed point iteration to this new mapping, the 
convergence speed is improved in a well-defined sense, and 
the computational complexity is not unduly increased because 
only one additional matrix-vector multiplication per iteration 
is required. As exemplary applications of the above results, 
we consider the problems of power and load estimation in 
OFDMA-based systems 0, E). 

This study is structured as follows. In Sect. |II] we review 
basic results in convex analysis and in interference calculus. 
The material in Sect. |II] can now be considered standard, but 
we also show a simple proof of the fact that positive concave 
functions are standard interference functions. In Sect. im we 
relate some results in 0 (used to compute lower bounds 
for load in LTF network planning) to standard results on 
recession functions in convex analysis. The relations are used 
in Sect. |IV] to derive conditions for the existence of fixed 


points of general positive concave mappings. We also pro¬ 
pose novel low-complexity iterative methods that improve the 
convergence speed of the standard fixed point algorithm. In 
Sect. |V]we revisit the problems of load and power estimation in 
OFDMA-based networks, and we show how the novel results 
and algorithms proposed here can be used in these concrete 
applications. 

II. Preliminaries 

In this study, we use the following standard definitions. By 
{x, y) for arbitrary x € K.^ and y € K.^, we denote the 
standard inner product {x,y) := x*y. Its induced norm is 
given by ||a;|| := ^yjx^x). The set I := A^} is the 

set of indices of the components of vectors in K.^. The oo- 
norm of a vector x = [xi,..., xat] is the norm given by 
ll^lloo := maxigx|a;i|. We define by the fcth standard 
basis vector of K.^. Vector inequalities should be understood 
as component-wise inequalities, and we define := [0, oo[^ 
and := ]0,cx)[^ (the superscript is omitted if A^ = 1). 
The set of positive integers is denoted by N := {1, 2,...}. 
The spectral radius of a matrix M G is given by 

p{M) := max{|Ai|,..., |Aiv|}, where Ai,...,AAr are the 
eigenvalues of the matrix M. The component of the Ah row 
and fcth column of a matrix M is denoted by [M]i^k- For a 
vector a: = [a:i,..., xn] G R^, the matrix diag(a;) G R^^'^ 
is a diagonal matrix with [diag(a;)]i j = Xi. 

Concave functions and standard interference functions play 
a crucial role in this study, so we review below basic definitions 
and known results that are extensively used in the next sections. 

Definition 1 (Convex set) A set C C R^ is said to be convex 

if 

(fJx G C)(yy G C')(Va G ]0,1[) ax + {1 — a)y G C. 

Definition 2 (Concave functions) We say that f : C —>■ R U 
{—oo} is a concave function if C C R^ is a convex set and 

(Va: G dom f){yy G dom f)(fJa G ]0,1[) 

f{ax + (1 - a)y) > af{x) + (1 - a)/(y), 

where dom f := {x G C \ f{x) > —oo} is the (effective) 
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domain of f. [[| 

Concave functions / : C —>■ R U {— cxd} with C C 
for which there exists at least one vector x € C satisfying 
f{x) > —oo are called proper concave functions. If for every 
sequence C C converging to an arbitrary vector x G 

C, we have limsup„_^g^ f{xn) < f{x), then we say that the 
function / is upper semicontinuous (on C). 

Every concave function / : C ^ RU {— cxd} can be related 
to a convex function by —/ : C —R U {00} (—/ takes the 
value CXD whenever / takes the value —cxd), so the following 
results on concave functions can be directly deduced from 
standard results on convex functions found in the literature 

im, Ida. 

Definition 3 (Superdifferentials and supergradients) Let 
f : C —> R U {—00} be a concave function with 

0 7 ^ C C R^. The superdifferential of f at x G dom / is the 
set given by 

df{x) ■= {m e R^ I (Vy gC) {y- x,u) + f{x) > f{y)} . 

If X ^ dom /, then we define df{x) := 0. A vector g G 
df{x) is called a supergradient of f at x. The domain of 
the superdifferential df is the set given by dom df := {x G 

R^ I df{x) 7 ^ 0}. 

In this study, if the point at which a supergradient is 
selected needs to be explicitly known, then we often use the 
notation g{x) G df{x) to denote an arbitrary choice of the 
supergradient at x. 

As a particular case of lfT3l Corollary 16.15], we have the 
following: 

Fact 1 Let f : —>■ R++ be concave. Then the superdiffer¬ 

ential df{x) is nonempty for every x G R^+. 

The proposed acceleration methods are based on the concept 
of recession functions (or asymptotic functions) in convex 
analysis. 


’in the literature, when a concave function / is allowed to take the value 
— 00 , assuming that C = is a common practice. If C is a proper subset 
of , we can define f{x) = —00 if cc € MJ^\C to extend / from C to 
R’'^. By doing so, the effective domain is preserved. However, for notational 
convenience later in the text, we do not necessarily adhere to this convention, 
and we allow C to be a strict subset of 


Definition 4 (Recession or asymptotic functions) Let f 
R^ —>■ R U {—cxd} be upper semicontinuous, proper, and 
concave. We define as its recession or asymptotic function at 
y G R^ the function given by (see U2\ Ch. 2.5] 4731 p. 152] 
for the standard definition for convex f): 

rw ^ j rt r t ^ r M “ fix) 

(Vai G dom/) /oo(y) := lim ---. 

h—¥oo ri 

(NOTE: The above limit is always well defined. We assume 
that it can take the value —cxd.) 

Fact 2 7//: R^ ^ R U {—cxd} is a proper, upper semicon¬ 
tinuous, and concave function, then for every y G dom/ we 
have M2\ Corollary 2.5.3] 

fooiy) = lim h/(h"^y), ( 1 ) 

/1-S-0+ 

and the above is valid for every y G R'^ if 0 G dom /. 

Fact 3 Let f : R^ —>■ R U { — 00 } be proper, upper semicon¬ 
tinuous, and concave. Then M2\ Proposition 6.5.1] 

/oo(y) = inf{(y,y) | a; G dom 5/, gG df(^x)}. 

Many estimation and optimization tasks in communication 
networks can often be posed as systems coupled by standard 
interference functions, which we define below. 

Definition 5 (Standard interference functions and mappings 
ffTj) A function f : R)^ —^ R+-l_ is said to be a standard 
interference function if the following properties hold: 

1) (Scalability) (ffx G R+) (Va > 1) af{x) > f{ax). 

2) (Monotonicity) (ffxi G R+) (Vx 2 G R+) Xi > X 2 ^ 
f(xi) > f(x 2 ). 

Given N standard interference functions fi : R++, 

i = 1,..., TV, we call the mapping T : R)j} : X I—>■ 

[/i(a;),..., f]\r(^x)] a standard interference mapping. 

Fact 4 (Properties of interference mappings m) Let T : 
R^ —>■ R+_|_ be a standard interference mapping. Then the 
following holds: 

Fact 4.1 Let Fix(T) := {x G R++ | T^(x) = x} be the set 
of fixed points of T, then Fix(T) is either an empty set or a 
singleton. 

Fact 4.2 Fix(T) 0 if and only if there exists x' G R'^ such 
that T{x') < x'. 
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Fact 4.3 If Fix(r) f 0, then it is the limit of the sequence 
generated by Xn+i = T{xn), where Xi G is ar- 
bitrary^If Xi satisfies T(xi) > Xi (resp. T(xi) < Xi), then 
the sequence is monotonically increasing (resp. monotonically 
decreasing) in each component. In particular, monotonically 
increasing sequences are produced with Xi =■ 0 . 

The focus of this study is on (positive) concave functions, 
which as shown below are a subclass of standard interference 
functions. 

Proposition 1 Concave functions f : M.^ -G- M++ are 
standard interference functions: 

Proof: We need to prove that concave functions / : 

M++ satisfy the scalability and monotonicity properties in 
Definition |3 

(Scalability) Let p > 1 and x G R+ be arbitrary. By 
concavity of /, for every a G ]0,1[, we have f{apLx) = 
f{afix+{l — a)0) > af (fix) + {1 — a) f {0). In particular, for 
a = 1 /p,, we conclude from the last inequality and positivity 
of / that 

fix) > -fipx) + ("l - /(O) > -finx), 

p \ P / P 

which proves the scalability property. 

(Monotonicity) Let {xi,X 2 ) G x satisfy X 2 > xi. 

As a result, Xi + p,{x 2 — xf) G for every p > 0. From 

the definition of concavity, we also have 

(Va G ]0, l[)(Vp > 0) 

/ ((1 - a)£Ci +a{xi+ p{x 2 - tCi))) > 

(1 - a)fixi) + af {xi + p(£C 2 - Xi)) . 

In particular, for an arbitrary p > 1 and for a = 1/p., we 
obtain from the positivity of / that 

f{x2) > f 1 - fixi) + -f {xi + p{x2 - Xi)) 

> f{xi) --fixi). 

The inequality /(a? 2 ) > /(®i) “ (1 /m)/(^i) is valid for 
every // > 1 , so we can take the limit as // goes to infinity to 

^In finite dimensional spaces, all norms are equivalent. Therefore, conver¬ 
gence of a sequence {aJnjneN C to a point E does not depend 
on the choice of the norm; i.e., limn-^oo ||®n — ®*||a = 0 for any norm 


conclude that 

fix 2 )> lim ff{xi) --f{xi)'] = f{xi). 

■ 

As every result stated in this section. Proposition [T] can 
be considered standard (see 0 and the references therein). 
Nevertheless, we have decided to include a simple proof of this 
proposition because similar statements can often be found in 
the literature without proof. Furthermore, some partial proofs 
available in the literature make implicit assumptions such as 
the existence of the supergradients on the boundary of the 
domain M.^ and/or the strict concavity of the functions. We 
emphasize that these assumptions are not required. As an 
example of a positive concave function (and hence a standard 
interference function) not satisfying these two assumptions, we 
have 

f 1 , if a; = 0, 

/ : M+ —^ M++ : a; !->■ < 

I 2 , otherwise. 

To characterize the existence of fixed points of affine stan¬ 
dard interference mappings, we can use the following fact: 

Fact 5 /E] Theorem A. 16] For an arbitrary matrix M G 

//■ piM) < 1, then converges and [I — 

M)-i = J + 

Fact 6 /E] Theorem A.51] Let M G be a non¬ 

negative matrix, and let p G R//+ be arbitrary. A sufficient 
and necessary condition for the system x = p Mx to have 
a (strictly) positive solution x G R++ is p(Vf) < 1. 

We end this section with a very simple statement that is 
used later to clarify an argument in Sect. IVl 

Remark 1 Let M G be arbitrary and D G R^^-^ be 

an invertible matrix. Then the eigenvalues of the matrices M 
and DMD~^ are the same (which in particular implies that 
p{M) = p(DMD^^)). 

Proof: Assume that a; is a right eigenvector associated 
with an eigenvalue A, and define y := Dx. As a consequence, 

Mx = Aa; 4=> MD~^y = XD~^^y DMD~^^y = Ay, 


and the result follows. 
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III. Component-WISE ineimum oe supergradients oe 
POSITIVE concave EUNCTIONS 

The main objective of this section is to propose two simple 
techniques for computing the component-wise infimum of 
supergradients of concave functions (c.f. Proposition |2] and 
Proposition O. These techniques are motivated by the fol¬ 
lowing application. In load estimation problems in wireless 
networks, the values taken by partial derivatives of functions 
related to the load coupling among base stations attain their 
infimum asymptotically as we move to infinity in the direction 
of a basis vector a, m- This observation has given rise to 
efficient techniques for the computation of lower bounds for 
the load in that very particular application domain na, and 
extending these results to a more general class of concave 
functions is highly desirable for other applications such as 
power estimation in networks. 

By using the concept of recession or asymptotic functions, 
we show below that the above-mentioned asymptotic result 
can be generalized to all positive concave functions, even if 
the functions are not differentiable, in which case we use 
supergradients instead of gradients (c.f. Proposition O. We 
can further show that the component-wise infimum taken by 
the supergradients can be easily obtained by means of simple 
schemes that do not require the computation of supergradients 
(c.f. Proposition Ell. These infimum values are used later by 
the proposed acceleration schemes to compute fixed points of 
positive concave mappings, and they can also be used to obtain 
a certificate that the mapping does not have a fixed point. We 
start by formalizing some simple properties of supergradients 
of concave functions. 

Lemma 1 Let f : —5- R++ be an upper semicontinuous 

concave function. Then the following holds: 

Lemma 1.1 All supergradients of f are non-negative vectors; 
i.e., 

(Va; G dom df) (Vg G df{x)) g > 0. 

Lemma 1.2 Let x G and k G I be arbitrary and assume 
that X -f /iCfe G dom df for every h > 0. Then 
(V/i > 0) (Vg' G df{x)) (Vg" G df{x -f hek)) 

0<9k< 9l 

where [gj,..., gJv]* — d' [g'{, ..., ■= a" ■ 


Lemma 1.3 As in Lemma 17.21 let x G and k G I be 
arbitrary and assume that x-\-hek G dom df for every h > 0. 
Then 


ifJh > 0)(Vg(a: -f hek) G df{x + heu)) 


9k H” — 7 5 (^) 


where [gi{x -f hek ),..., g7v(ic + hek)Y ■= g{x + hek). 


Proof: 

1) We prove the result by contradiction. Assume that there 
exists a supergradient g =: [gi,...,g^v]* G df{x') at 
some point x' G dom df such that g* < 0 for an 
arbitrary i G {l,...,iV}. We know from the definition 
of supergradients that 

(Vg G R+) f{y) < fix') -f g\y - x'). 

In particular, for tt : M —>■ : h x' hei, we 

obtain 

fiu{h)) < fix') -b g'iuih) - x') = fix') -b g^h. 

Now, since g^ < 0 by assumption, we obtain /(tt(li)) < 
0 for an arbitrary h > /(a;')/|gi|, which contradicts the 
positivity of the range of the function / : -G M++. 

This proves Lemma o 

2) By Definition 13 for arbitrary a;i,a ;2 G dom df, we 
have fixi) < fixf) -b glixi - Xq) and fix^) < 
fixi) -b g\ixo - xi), where g^ G df ixo) and g^ G 
dfixi) are arbitrary supergradients. Summing these 
two inequalities yields 

iai-9oY ixi-xo) <0. (3) 

In particular, for Xi = x hek and Xq = x, we 
have 0 < a;i — a;o = bek ^ 0 , and we can set 
go = g' G dfix) and g^ = g" G dfix -b hek). Using 
these particular choices for Xg, Xi, g^, and g^ in 0 , we 
obtain g'f < g'^. Non-negativity of g'f has been proved 
in the first part of the lemma. 

3) Use a:' = a; in the supergradient inequality fix') < 
fix -b hek) + aix -b hekfix' - x- hek). 


We can now show an efficient scheme to compute the 
element-wise infimum of supergradients. 
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Proposition 2 Let S := ^/(®) ®/ 

supergradients of an upper semicontinuous concave function 
f : ®^++- For each fc € X, define 

9k ■= inf {ff/c G R I g = [ 51 ,..., ffiv] G S'} € R+, (4) 

then we have 


(VfcGJ) gl= lim hf{h ^e^). 

h-^0+ 

Proof: 

We have g]! > 0 as a direct consequence of Lemma fTTI Now 
consider the standard extension / : R^ — 1 ^++ U {—00} of 
/ : R^ —>• R++ given by 

f fix) if a; G R? 

/(^) = r^ + 

I —00 otherwise, 

and let fc G X be arbitrary. By construction, / is upper 
semicontinuous, proper, and concave. Furthermore, dom df = 
dom df, dom / = dom/, and df{x) = df{x) for every 
X G dom df. By Fact |2 we have = fooi^k) (see 
Definition |4]i, and the result now follows from Fact |2] ■ 

Next, we show an alternative means of computing g^ in 
(Si. This alternative method has been used in m for a very 
particular concave function appearing in load estimation in 
LTE networks (see Sect. IV-Ab . 

Proposition 3 Let x G and k € I be arbitrary. In addi¬ 
tion, assume that f : R:^^ —>■ R++ is an upper semicontinuous 
concave function and that x-\-hek G dom df for every h > 0. 
Define by [gi{x + hek),. . • ,5Ar(a; + hek)] := gix + hsk) G 
dfix + hek) arbitrary supergradient at x -\- hek. Then 

lim gkix + hek)= g*k> 0, (5) 

h—foo 

where gf. is defined in (0. 

Proof: Let /c G X be arbitrary. It follows from Lemma [L2] 
that, irrespective of the criterion we use to select a super¬ 
gradient gix + hek) G dfix + hek), its fcth component 
gkix + hek) should be monotonically non-increasing as h 
increases (and lower bounded by 0). As a result, the limit 
lim^_,.oo gkix + hek) exists. By definition, is the infimum 
of the A:th component of all supergradients, hence we have that 

gl < lim gkix hek) 

h—¥oo 


for any choice of x and k satisfying the assumptions of the 
lemma. Using (|2l) in Lemma [T3] and the definition of recession 
functions, we deduce 


< lim gkix -f hek) 

h—¥oo 


< 


lim 

h—¥oo 


fjx-i-hek)- fjx) 
h 


fooi^k)- 


The result now follows by noticing that fooi^-k) = gt by 
Fact [ 3 ] (Non-negativity of g^ is immediate from Lemma ll.lh 


IV. Acceleration algorithms eor positive 

CONCAVE MAPPINGS 

Having two efficient methods to compute the component¬ 
wise infimum of supergradients of concave functions, we can 
now proceed with the study of general concave mappings. 
To avoid unnecessary technical digressions, we do not deal 
with concave functions / : R^^ ^ R++ that are not upper 
semicontinuous. To formalize this assumption, we use the 
following definition; 

Definition 6 (Positive concave mappings) We say that T : 
M.^ —>■ R++ is a positive concave mapping if it is given by 

Tix) ■.= [fiix),...,fNix)W ( 6 ) 

where all functions fi : R:^^ —>■ R++,... , /at : R+ —>■ R++ 
are concave and upper semicontinuous. 

By Proposition [T] we know that positive concave mappings 
are standard interference mappings. The remaining of this 
section has the objective of investigating the following prob¬ 
lems associated with a positive concave mapping T (which, as 
shown in Sect. [V] are problems that need to be addressed in 
many network planning and optimization tasks): 

PI) Verify whether T has a fixed point by using computa¬ 
tionally efficient algorithms. 

P2) Improve the convergence speed of the standard iteration 
in Fact 14.31 to obtain the fixed point of T (if it exists). 

A. Conditions for the existence of fixed points of positive 
concave mappings 

To address problem PI), we use the concept of lower 
bounding matrices, which we define as follows: 





7 


Definition 7 The lower bounding matrix of a positive concave 
mapping T : ^ R^_|_ : x i-^ [fi{x),..., f]\f{x)Y is the 

non-negative matrix M £ with its ith row and kth 

column given by 

[M]i^k ■= inf {gfe e K I [gi,... ,g7v] e 54 £ IR+, (7) 
where Si := UcceR" 9fi{x). 

Note that Proposition |2 and Proposition [3] show two simple 
techniques to compute each component of lower bounding 
matrices. 

Example 1 (Construction of lower bounding matrices with 
the results in Propositions \2\ and |5]l Let the functions fi : 
R^ —t IR++,..., /at : —>■ M++ be concave and upper 

semicontinuous. Using Proposition |2] we can compute the 
lower bounding matrix M of the mapping T : ^ R^^+ : 

X [fi{x),.. .jNix)f by 

M = 

lim;j^o+ hfi{h-^ei) ■ ■ ■ lim^^o+ hfi{h-^eN) 
lim;,_5 .o+••• lim,,^o+ 

( 8 ) 

Equivalently, by fixing x' £ R^+ arbitrarily, we can also 
compute the lower bounding matrix M with the results in 
Proposition |5] and Fact Q] as follows: 

M = 

lim/j^oo gKx' + hei) ■ ■ ■ limA-j-oo glfix' + Hen) 

! *. I ? 

gi (x'hei) ■■■ limh^oo g^ix'heN)_ 

(9) 

where we denote by g],{x) the kth element of a supergradient 
of ft at X e R^+ (i.e., [gl{x),.. .,g}^{x)Y £ dffix)). 

Proposition |2] and Proposition [3 also show that the lower 
bounding matrix is non-negative. The name “lower bounding 
matrix” stems from the fact that this matrix is constructed 
with component-wise lower bounds of supergradients. Lower 
bounding matrices can also be used to construct affine map¬ 
pings that serve as lower bounds of their corresponding posi¬ 
tive concave mappings, in the following sense: 


Lemma 2 Let M be the lower bounding matrix of a positive 
concave mapping T : R^ ^ R^+ in accordance with 
Definition 0 Then 

(Vy £ R^)(Va; £ R^) x >y ^ T{x) > T{y) M{x — y). 

( 10 ) 

Proof: We prove the inequality for an arbitrary component 
of the mapping T; i.e., for the function fi, where z £ I is 
arbitrary. Let y £ R:!J( and a; > y be arbitrary vectors, and 
construct the sequence {xn := (l/n)l -t- x}n gn C R++- By 
Fact [T] we have a;„ £ dom dfi for every n £ N. From the 
definition of supergradients, we know that, for every n £ N, 

My)+gUxn-y) < Mxn), (ii) 

where y„ £ dffixn) is an arbitrary supergradient. By Defini¬ 
tion I 2 ] the zth row of M, denoted by rrii > 0 as a column 
vector, is the component-wise infimum of all supergradients of 
the function fi, hence 0 < < y„ for every n £ N. Using 

this last relation together with a;„ > y in (fTTT i. we deduce: 

(Vn £ N) 

fi{y) -f ml{xn -y)< fiiy) -t- y^(a;« - y) < fi{xn). 

By construction, lim„_>oo a;„ = a:. As a result, we con¬ 
clude from the continuity of affine functions and upper semi¬ 
continuity of fi that 

fi{y) + ml{x - y) = limsup(/,(y) -b m-(at„ - y)) 

n—foo 

< limsup/i(at„) < fi{x). 

n—foo 

■ 

The next proposition addresses problem PI) stated in the 
beginning of this section: 

Proposition 4 Let M be the lower bounding matrix of a 
positive concave mapping T : R;?^ —>■ R+_|_. A necessary 
condition for Fix(T) f 0 is p{M) < 1. 

Proof: Use y = 0 in (fTOl) to verify that the affine mapping 
Tl : R+ —R++ : x 1 —>■ T{0) Mx satisfies Tlx) > 
Ti^{x) > 0 for every x £ R:^^. Being an affine mapping, 
Tl is a positive concave mapping, hence it is also a standard 
interference mapping by Proposition [T] Now let x* £R?+ be 
the fixed point of the mapping T. By Lemma |2] we obtain: 


Tl[x*) < T{x*) = xfi 


( 12 ) 






which implies the existence of the (unique) hxed point of the 
mapping Ti^ by Fact|4]2. In other words, there exists a unique 
positive vector x G satisfying x = T(0) + Mx, and we 
know by Fact|6]that there exists a positive vector satisfying this 
last equality if and only if p{M) < 1 (recall that r(0) > 0 
and that M e by construction). ■ 

An immediate consequence of Fact |5] and Proposition |4] is 
the following useful result: 

Corollary 1 Let T : ®++ ^ positive concave 

mapping with Fix(r) ^ 0, and denote by M € its 

existing lower bounding matrix given by Definition 0 Then 
(I — exists, and it is a non-negative matrix. 

Proposition 0] is interesting in its own right because it 
enables us to certify that a given positive concave mapping 
has no hxed point. We only need to show that the spectral 
radius of its lower bounding matrix has spectral radius greater 
than or equal to one. This result is highly relevant in network 
optimization and planning problems. As already mentioned 
in the introduction, in these applications, the feasibility of a 
network design follows from the existence of the hxed point 
of a mapping that is constructed based on antenna tilts, power 
allocations, the position of base stations, etc. Optimization of 
the network performance (e.g., in terms of energy efficiency, 
capacity, coverage, etc.) over the joint set of all control 
parameters is typically an NP-hard problem. As a result, many 
optimization algorithms proposed in the literature are greedy 
heuristics that need a fast feasibility check of multiple network 
conhgurations at each iteration a. Proposition 0] opens up 
the door to the development of efficient and fast methods 
for excluding many infeasible network conhgurations from 
consideration, which can signihcantly accelerate the overall 
optimization process. 

We emphasize that the converse of Proposition 0] does not 
hold in general. There are mappings for which the lower 
bounding matrix has spectral radius strictly less than one, and 
yet mappings do not have a hxed point (see the application 
in Sect. IV-Bb . Therefore, to characterize the existence of a 
hxed point based on the spectral radius of the lower bounding 
matrix, we need additional assumptions on the mapping. The 
next proposition shows a particularly useful assumption that is 
satished in load estimation problems (see Sect. IV-AI and ||6l 
for a particular application of this proposition). 


Proposition 5 Let T : —>■ K.++ be a positive concave 

mapping with lower bounding matrix M satisfying p{M) < 1. 
In addition, assume that 

(3y GM^+)(Va; T{x) <y + Mx. (13) 

Then the mapping T has a fixed point. 

Proof: Let y' G K^+ be a vector satisfying T{x) < 
y' + Mx for every x G R;!^. By Fact. |6] we know that x' := 
{I — M)~^y' is a strictly positive vector. Therefore, 

T(x') <y' + Mx' = (I - M)x' + Mx' = x', 

and the above implies that Fix(T) ^ 0 by Fact. 14.21 

■ 

B. Acceleration techniques for positive concave mappings 

We now turn our attention to problem P2). To address this 
problem, we use the concept of accelerated mappings, which 
we dehne as follows: 

Dehnition 8 (Accelerated mappings) Let T : —>■ M+_|_ be 

a positive concave mapping and M be its lower bounding 
matrix. If p{M) < 1, the accelerated mapping Ta : —>■ 

of T is the mapping given by: 

TaIx) := (I - M)-\T{x) - Mx). (14) 

To see that the codomain of Ta in the above dehnition is 
indeed note that, by Lemma |2l we have that T{x) — 

Mx > T(0) > 0 for X G R+. Now use Fact |5] to conclude 
that (J — M)~^ [T{x) — Mx) is a (strictly) positive vector 
for every x G R;!^. 

The next lemma shows an alternative way to compute 
accelerated mappings. The main advantage of this alternative 
expression is computational. We have to perform only one 
matrix-vector multiplication. 

Lemma 3 Let Ta : R^ ^ R(^+ be the accelerated mapping 
of the concave mapping T : R^ ^ R^+, where we assume 
that the lower bounding matrix M of T satisfies p{M) < 1. 
Then Ta in l\14[ can be equivalently expressed as 

Ta{x) = (J - M)-\T{x) -x)+x. 


(15) 
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Proof: Recalling that the matrix I — M is invertible as a 
direct consequence of Fact|3 we deduce: 

T^ix) = {I- M)-\Tix) -x)+x 

■^{I - M)Ta{x) = T{x) -x + {I- M)x 
O (I - M)Ta{x) = T{x) - Mx 

O Ta{x) = {I- M)-\T{x) - Mx). 

■ 

Positive concave mappings and their corresponding ac¬ 
celerated mappings have many common characteristics. In 
particular, they are both standard interference mappings, and 
they have the same fixed point, as shown below. 

Lemma 4 Assume that T : > R++ is a positive concave 

mapping with lower bounding matrix M satisfying p{M) < 
1. Then the accelerated mapping Ta ■ R+ —>■ R++ of T 
is a standard interference mapping. Furthermore, Fix(r) = 
Fix(rA). 

Proof: By p{M) < 1, the matrix inverse (J — M)~^ 
exists, and it is a non-negative matrix (Fact |5]). As a result, 
each component of the mapping T'{x) := (J — M)~^T{x) 
{x G is a positive sum of concave functions, hence 
the resulting function is also concave. Observing that linear 
functions are both concave and convex, we verify that each 
component of Ta{x) = T'{x) — [I — M)~^x -f a; > 0 is 
a positive sum of concave functions, hence Ta is a positive 
concave mapping. Proposition [T] now shows that Ta is a 
standard interference mapping. Consequently, Ta has a unique 
fixed point, if it exists (Fact 14.1b . If x* G Fix(r), then 
Ta{x*) = {I - M)-\T{x*) - Mx*) = (/ - M)-^{x* - 
Mx*) = (J — M)~^{I — M)x* = X*. The converse is 
also immediate. Note that T(x) = {I — M)Ta{x) + Mx, 
hence T(x*) = x* if Ta{x*) = x*, and we conclude that 
Fix(T) = Fix(TA). ■ 

The practical implication of Lemma|4]is that, to compute the 
fixed point of a positive concave mapping T, we can instead 
compute the fixed point of its accelerated version Ta by using 
the standard iteration Xn+i = TA{xn) shown in Fact |4] In 
many applications, having a monotone sequence is 

desirable, and a sequence of this type can be constructed with 
the standard fixed point iteration by starting the iterations from 
a;i = 0 (see Fact 14.3b . For example, in network planning and 
optimization tasks, the fixed points of the concave mappings 


are estimates of the power allocation or of the load at the base 
stations ii-ii. Therefore, even if the mapping has a fixed 
point, the network design is invalid if the power or load of any 
base station exceeds its physical limit. If the iterative algorithm 
produces a monotonically increasing sequence, we obtain a 
certificate that the design is invalid as soon as any element of 
the vector sequence exceeds its limit. It is particularly in these 
cases that the standard iteration with the accelerated mapping 
Ta converges faster than the standard iteration with the original 
mapping T, in the following sense: 

Definition 9 (Faster convergence) Let {a;„}„gp^ C and 
{yj C M.^ be two sequences converging to the same 
vector u* G We say that the sequence C 

converges faster than {y„}neN C if\\xn-u*\\ac, < ||y„- 
■*^*1100 for every n G N. 

With the above definition, we can now formally state the 
improvement obtained by using Ta instead of T with the 
standard iteration in Fact 14.31 

Proposition 6 Assume that T : > R+a is a positive 

concave mapping with lower bounding matrix M G R^^'^ 
satisfying p(M) < 1. Let Ta : R^ —>■ R++ be the accel¬ 
erated mapping of T. Consider the following two sequences: 

:= r(a;^)}„gN and := TA(a;")}„gN- Assume 

that both sequences start from the same vector u G R^/ 
i.e., u = x'l = x”. If is monotonically increasing 

(resp. monotonically decreasing) in each component, then the 
following holds: 

Proposition 6.1 is monotonically increasing (resp. 

monotonically decreasing) in each component. 

Proposition 6.2 a;" > ( resp. a;" < at^ j for every n G N. 

Proposition 6.3 If the mapping T has a fixed point ( which, 
in particular, it is automatically guaranteed if {a;^}„gN A 
monotonically decreasing in each component), then {a:"}„gN 
converges faster than {a;(j}„gN to x* G Fix(T), in the sense 
of Definition [9] 

Proof: We prove the proposition only for monotonically 
increasing sequences (in each component). The proof for 
monotonically decreasing sequences can be obtained in a 
similar fashion by reversing all inequalities. 
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1) Recall that, by Lemma E] Ta is a standard interference 
mapping, so, in light of Fact 14.31 we only need to prove 
that Ta{u) > It if T{u) > u. 

By assumption, T{u) — tt > 0 and p{M) < 1. In 
particular, by using Fact 0 and non-negativity of M, 
the last inequality implies that (7 — M)~^ is a non¬ 
negative matrix. Consequently, from (fTSl l. we deduce; 

Ta{u) = (7 — M)~^(T{u) — u) + u > u. 

2) We show the result by using induction. Assume that 
tc" > x'^ for a given n S N. From the definition of the 
mapping Ta in (fT4l) . we deduce: 

Ta{x^) = = T{x^) -f — x^). (16) 

We have already proved that {a;"}„gN is monotonically 
increasing with the assumptions of the proposition 
(hence - < > 0), M is a non-negative ma¬ 

trix, and T is a mapping satisfying the monotonicity 
property of standard interference functions. Using these 
observations in (fThl) . we verify that; 

^n+i ~ T^Aix^) > T(x^) > T{x^) = 

The above arguments are valid, in particular, for n = 1, 
because x'l = x'^= u by assumption. 

3) First recall that both {£c^}„gN and {a;"}„gN converge 

to the uniquely existing fixed point x* £ Fix(T) 
(Lemma |4] and Fact 14.31 ). The desired result ||a;" — 
a;* I loo < — ic*||oo, valid for every n £ N, follows 

directly from Proposition |6l2. 

■ 

As an immediate consequence of Proposition |6] and Fact 14.31 
we have the following. 

Corollary 2 Assume that T : -A K.a_|_ is a positive 

concave mapping with x* £ Fix(T) ^ 0, and denote by 
Ta : —>■ R++ its corresponding accelerated mapping. 

Then the sequence {TjJ(0)}„gN converges faster to x* than 
the sequence {T”(0)}„gN, in the sense of Definition^ (we 
assume that both sequences start from the vector 0). 

Remark 2 Following Yates’ arguments ^7^ to prove the con¬ 
vergence of the iteration in Fact. 14.51 we can also argue 
that the proposed accelerated scheme is expected to be fast 
when the initial point is arbitrary. More precisely, assume 


that T : —> R+a is a positive concave mapping with 

a fixed point denoted by x* £ Ra+- Since this fixed point is 
strictly positive, for an arbitrary vector x £ R:^^ there always 
exists a > \ satisfying x < ax*. From Definition |5] we verify 
that T{x) < T{ax*) < aT{x*) = ax*. These inequalities 
imply that (see also Fact. 14.51) i) T"(0) < T^{x) < T"(aa:*) 
for every n £ N, ii) the sequence {r"(aa;*)}„gis} is mono¬ 
tonically decreasing, and Hi) {T"(0)}„gN is monotonically 
increasing. In other words, each term of the monotone se¬ 
quences {T”(0)}„gis} and {r"(aa;*)}„gN are, respectively, 
(element-wise) lower and upper bounds for each term of 
the sequence {T"(a;)}„gN. All the above arguments are also 
valid if we exchange T by its corresponding accelerated 
mapping Ta, and we note that lower and upper bounding 
sequences {TjJ(0)}„gN and {fr'^{ax*)}n&i for the sequence 
{T^(a;)}„gN converge faster to the fixed point x* when com¬ 
pared to the lower and upper bounding sequences {T”(0)}„gi!^ 
and {T”(aa;*)}„gN/or the sequence {T"(a;)}„gN. In other 
words, the sequence produced by a;^_|_]^ = Ta{x'.^) with 
x[ = u G R^ arbitrary has sharper element-wise bounds 
than the sequence produced by a:"_|_]^ = T(x'f) for the same 
starting point x'l = u. 

Remark 3 The price we pay to use the accelerated iteration 
x'^_^i = Ta{x'^) instead of using Xn+i = T{xn) is the need 
for a matrix-vector multiplication, if Ta is evaluated by using 
(Ell (assuming that the lower bounding matrix is not the zero 
matrix). Furthermore, a matrix inversion is required (or, for 
increased numerical stability, a matrix decomposition), but this 
operation needs to be done only once. One situation where the 
proposed scheme is particularly useful is when the evaluation 
of the mapping T is time consuming when compared to the 
matrix-vector multiplication. In this situation, for all practical 
purposes, the time to compute x'^ or a;„ is roughly equivalent 
for a given n £ N sufficiently small. However, for every n £ N, 
x'^ is guaranteed to be a better approximation of the fixed 
point of the mapping T than x„. This situation is common in 
network planning. 

V. Applications in Network Planning and 
Optimization 

We now apply the general results in the previous sections 
to two concrete estimation problems in LTE networks. First, 
we consider the load estimation task discussed in ii-ilil, 
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among other studies. Briefly, the objective is to determine the 
bandwidth required to satisfy the data rate demand of all users 
in the network, by assuming that the transmit power of all base 
stations is given. 

The second application we consider is the reverse of load 
estimation. For a given load allocation at the base stations, 
the objective is to estimate the power allocation inducing that 
load. This reverse problem has been motivated by the study in 
il, which has proved that using all available bandwidth is ad¬ 
vantageous from various perspectives, and, in particular, from 
the perspective of transmit energy savings and interference 
reduction. That study also proves that there exists a standard 
interference mapping having as its fixed point the solution of 
the power estimation problem, and the study in ||9l has shown 
that the interference mapping can take the form of a positive 
concave mapping. 


nonlinear equations 141-171: 




Vm = fM{t^,P), 


where 


(17) 


/, : Rfx ^ R++ 




Note that, for each fixed p and i £ A4, the function 

hp.i : —?> R++ : ix fi{ix,p) is concave, hence the 

solution of (fm with fixed p can be obtained by computing 
the fixed point of the positive concave mapping given by 0, 

Q 


Tp{v) := [hp,i{v ),..., (18) 

A. Load estimation 


We consider an LTE network with M base stations and N 
users represented by elements of the sets M. = {1 ,..., M} 
and M = {1,..., N}, respectively. The set of users connected 
to base station i S A4 is denoted by Mi, and the data rate 
requirement of user j S AA is given by dj > 0. The propagation 
loss between user j £ M and base station i £ Ml is denoted by 
gij > 0. Each base station i £ Ml has K resource units that 
can be assigned to users, and the transmit power per resource 
unit for each base station i £ Ml is pi >0. The reliable 
downlink data rate for each resource unit connecting base 
station i £ Mi to user j £ M is approximated by the following 
well-established interference-coupling model 0-0, EO): 




Blog2 



Pi9i,j 

EfeGAr\{j} ^kPk9k,j 



where is the noise power per resource unit, p = 
\pi,... ,Pm]* is the downlink power vector per resource unit, 
h> = [vi,..., i^mY is the load vector, and B is the bandwidth 
per resource unit. Here, the load Vi is fraction of the number 
of resource units in the time-frequency grid that users in the 
set Mi require from base station i. Eor fixed power allocation 
p £ KE_, the load is the solution to the following system of 


Therefore, all the theory developed in the previous sections ap¬ 
plies to this problem, and, in particular, the novel acceleration 
schemes for the computation of fixed points. Before proceeding 
with numerical examples of the acceleration schemes, we 
revisit known results related to this problem, and we show how 
the application-agnostic approaches developed in Sect. [Ill] and 
in Sect. |IV]can be used to reach these known results in a more 
convenient way. 

In particular, the authors of 0 construct a matrix by com¬ 
puting the values that the partial derivatives of the functions 
/ip.i,..., hp^M attain when a given component of the argument 
1 / of these functions goes to infinity. It has been shown in ifTbl 
that the system of nonlinear equations in (fTTl i has a solution if 
and only if the spectral radius of this matrix proposed in 0 
is strictly less than one. Using the terminology and results in 
Sect. Hn] and in Sect. IIVI we note that the matrix suggested 
in 0 is a particular case of a lower bounding matrix in 
Definition 0 constructed with the technique in Proposition 0 
The fact that the spectral radius of this lower bounding matrix 
gives sufficient and necessary conditions to characterize the 
existence of a solution of the nonlinear system is a direct con¬ 
sequence of the application-agnostic results in Proposition |4] 
and Proposition 0 

To be more precise, we can use (|9]l to construct the lower 




12 


bounding matrix Mp of the mapping Tp as follows: 


Mp = 


limh^oo gl{t^'+ hex) ■ 

■ ■ lim,i^oo ffM 

gf^+ hex) ■ 

• • limzi^oo 5 m + hsM ) 


where v' e is arbitrary and is the fcth component 

of a supergradient of the function hp^i at an arbitrary point 

u = [z2 i,..., z/m]*- By noticing that the function hp^i is 

differentiable in the interior of its domain, glii') is simply 

d 

the partial derivative —hp Av) for every i> S M.¥,. As a 
result, we can verify that the lower bounding matrix of the 
mapping Tp is given by Mp = diag(p)“^iVf'diag(p), where 


fo, 

[M']i^k = S ^ ln{2)djgkj 


if i = k 
otherwise. 


(19) 


(We can also obtain (fl^ by constructing the lower bounding 
matrix with the approach in ®.) 

By RemarklU p{Mp) < 1 is equivalent to p{M') < 1, and 
we note that M' does not depend on the power allocation p, 
a fact originally stated in ifTbll . Therefore, to verify whether 
the mapping Tp has a fixed point by using the results in 
Proposition |4] and Proposition |5l we can compute p{M') 
instead of p{Mp). In other words, knowledge of p{M') 
is sufficient to determine whether the system of nonlinear 
equations in ( fTTl i has a solution, as already stated in IfT^ for 
this particular application. 

Having the lower bounding matrix in closed form, we 
can now proceed to the numerical evaluations of the novel 
acceleration schemes. In the simulations we show here, we 
compare the accuracy of the load estimates generated by the 
standard iteration Un+i = Tp{un) with its accelerated version 
= Tpp^{v'^). Table IJ lists the main parameters of the 
network. 

The figure of merit used in the comparisons is the expected 
normalized mean error (NME), which we define by 


eNME(i^) := E[\\u - zz*||/||iz*||], (20) 

where v* G Fix(Tp). We approximate the expectation operator 
by averaging the results of 100 runs of the simulation, and 
in each simulation the positions of the users (and hence the 
propagation loss) are the random variables. All iterations start 
from the zero vector, and networks where the corresponding 


TABLE I. Network parameters of the simulation 


Parameter 

Value 

Caiiier frequency 

900 MHz 

Number of resource units (K) 

25 

Transmit power per resource unit {pi, Vi G A4) 

1.6W 

System bandwidth {K ■ B) 

5 MHz 

Noise power spectral density 

-145.1 dBm/Hz 

Propagation model 

Okumura-Hata, urban 

Antenna height of base stations 

30m 

Antenna height of the users 

1.5m 

Number of users (A^) 

200 

Number of base stations (M) 

25 

Data rate of each user {dj , \/j G A/”) 

768 kbps 

Dimension of the field 

2500mx 2500m 

User distribution 

Uniformly distributed at random 

Base station distribution 

Uniformly distributed 



Fig. 1. NME of the load estimate as a function of the number of iterations. 
Confidence intervals (95%) have been computed, but they are not visible in 
the figure. 


concave mapping does not have a fixed point are discarded. 
Therefore, the expectation in (|20] | is conditioned to the fact 
that spectral radius of the lower bounding matrix is strictly 
smaller than one. 

Fig.m shows results obtained by using the iterative scheme 
in Fact l4.3l with the original mapping Tp and with its proposed 
accelerated version Tp^. We verify that the mapping Tp^ 
requires fewer iterations than Tp to obtain a given numerical 
precision, which is an expected result by considering Proposi¬ 
tion |6] 


B. Power estimation 

We now turn our attention to the problem of power esti¬ 
mation in LTE networks. The objective is to solve (dill for 
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p = [pi, • ■ ■ ,PmY '''ith the load v € being the fixed 

parameter. It is shown in El that the solution of this nonlinear 
system is the fixed point of the positive concave mapping given 
by T^{p) ■= [hu.i{p), ■■■, h^,M(.P)Y’ where 


(P^ 


dj 


KAp) ■= ' 


— 1? -7-T’ if P* 7^ 0 

dj In 2 / ^ 2 \ 

-I/. ^ J : 


otherwise. 


By using ® to construct the lower bounding matrix 
of the mapping T^, we deduce: 


= 



■■ lim^^o+ 

lim3,^o+ xh^^M{x~'^ei) ■ 

• ■ lim3;_>o+ xK,M{x~^eM) 


= diag(i^) ^M'diag(i>'), 


where Al' is the same matrix defined in (fl^ . (The same result 
can be obtained by using Proposition [3] to construct the lower 
bounding matrix, but here applying Proposition|2]is easier than 
applying Proposition |3]) 

From Proposition 0] and the definition of Al^y, we conclude 
that a necessary condition for existence of the fixed point of 
Tiy is p{M') < 1, which is the same requirement for the 
existence of the fixed point of Tp in ( fTSl ). However, there 
is a fundamental difference between these two mappings. As 
proved in ifTfil . p{M') < 1 (note: this spectral radius does not 
depend on v) is both a sufficient and necessary condition for 
the existence of the fixed point of Tp. In contrast, the study 
in fSl has shown that the existence of the fixed point of T^, 
also depends on v. Therefore, T^, is an example of a mapping 
proving that the converse of Proposition 0] does not hold in 
general. 

We now turn the attention to the acceleration schemes in this 
particular application. We use the same network considered 
in the load estimation task. The desired load i/ is obtained 
by solving (fTTI) with the power fixed to the value shown in 
Table |I] Then we solve the reverse problem; we compute 
the power shown in Table U by using the standard iteration 
Pn+i = Tu{Pn) ^nd its accelerated version p'^+i = Ti,j^{p'.^). 
Both algorithms start from the zero vector. The normalized 
mean error is again used as the figure of merit (which in 



Fig. 2. NME of the power estimate as a function of the number of iterations. 
Confidence intervals (95%) have been computed, but they ai'e not visible in 
the figure. 


this application is defined by replacing the load vector by 
the power vector in (l20l) ). We can see in Fig. |2] that in 
this application the proposed acceleration scheme once again 
provides us with clear advantages over the standard iterative 
approach, in accordance to the analysis in Sect. IIV-BI 

VI. Conclusions 

We have shown that the results in m for the construction 
of lower bounding matrices in a very particular application 
domain can be generalized to a large class of positive concave 
mappings where even differentiability is not required. More 
specifically, we proved that positive concave mappings with 
nonempty fixed point set can be associated with a non-negative 
lower bounding matrix having spectral radius strictly smaller 
than one. By imposing additional assumptions on the mapping, 
having spectral radius strictly smaller than one also implies 
the existence of the fixed point of the concave mapping. 
We also demonstrated that the lower bounding matrix can 
be constructed with two simple and equivalent methods, and 
this matrix can be combined with its generating concave 
mapping to build a new mapping that preserves the fixed 
point. The standard fixed point iterations applied to this new 
mapping typically requires fewer evaluations of the original 
mapping to obtain an estimate of the fixed point for any given 
precision. The additional computational complexity of this 
novel approach is very modest. In the tasks of load and power 
estimation in LTE networks, where we are mostly interested 
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in the precision of the estimates after a limited number of 
iterations, numerical examples show that the improvement in 
convergence speed obtained with the proposed method can be 
substantial. 
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